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Abstract 

Quantum phases and phase transitions of weakly- to strongly-interacting 
bosonic atoms in deep to shallow optical lattices are described by a sin- 
gle multi- orbital mean-field approach in real space. For weakly- inter acting 
bosons in ID, the critical value of the superfluid to Mott insulator (MI) tran- 
sition found is in excellent agreement with many-body treatments of the Bose- 
Hubbard model. For strongly-interacting bosons, (i) additional MI phases 
appear, for which two (or more) atoms residing in each site undergo a Tonks- 
Girardeau-like transition and localize and (ii) on-site excitation becomes the 

excitation lowest in energy. Experimental implications are discussed. 
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Loading and manipulating cold bosonic atoms in optical lattices is a fascinating, rapidly 
growing branch of cold-atom physics, see, e.g., Refs. [1-3]. In a pioneering experiment, 
Griener et al. [1] have demonstrated the quantum phase transition from the superfluid (SF) to 
Mott insulator (MI) phase of a cold gas of ^^Rb atoms trapped by a three-dimensional simple- 
cubic-type optical lattice. More recently, the SF to MI transition has been demonstrated in 
an effective one-dimensional (ID) optical lattice [2]. In the SF phase, which has no excitation 
gap, atoms are free to move throughout the lattice and are associated with a coherent state 
of matter. The MI phase amounts for commensurate filling of the optical lattice and has 
an excitation gap associated with moving an atom from one site to an occupied neighboring 
site. The SF to MI transition (in deep optical lattices) involves weakly-interacting bosons 
and is well described by the Bose-Hubbard model [4-6] which assumes all bosons to occupy 
the lowest band of the lattice. Recently, the so-called Tonks-Girardeau gas, i.e., the strongly- 
interacting regime, was realized in an optical lattice [3] . Generally, as the interaction between 
atoms increases - in deep as well as in shallow optical lattices - additional possibilities open 
up for the trapped cold atoms which can now occupy higher bands. 

Our purpose in this letter is to explore quantum phases and excitations of cold bosonic 
atoms in optical lattices not accounted for so far. We will present an approach which is able 
to treat weakly- to strongly-interacting bosons in the entire range of deep to shallow optical 
lattices. To peruse the above, we would like to describe cold bosonic atoms in optical lattices 
directly in real-space., i.e., to provide their spatial wavefunction. Specifically, we will adopt 
our recently introduced multi- orbital mean-field [7,8] to cold bosonic atoms in optical lattices. 
As we shall see below, already for the standard case of weakly-interacting atoms we found 
a promising and intriguing result concerning the critical value of the SF to MI transition 
in ID. Although our approach is mean-field, it finds this value to be 3.855(7), in excellent 
agreement with density matrix renormalization group and other many-body calculations, see 
Ref. [9] and references therein. 

Our starting point is the many-body Hamiltonian describing N bosons in an optical 
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lattice, H — Y^f^i T{ri) -\-V{ri) + YH>j=iU{ri — Vj). Here, is the coordinate of the 
i-th particle, T(rj) and V{ri) stand for the kinetic energy and optical lattice potential, 
respectively, and U{vi — Vj) describes the pairwise contact interaction between the i-th and 
j-th atoms. 

As mentioned above, we are going to obtain a real-space, wavefunction picture of the 
quantum state of cold atoms in the optical lattice. How are we going to achieve that? To this 
end, we attach an orbital to each of the N atoms. The simplest choice is the Gross-Pitaevskii 
approach, for which all bosons reside in the same orbital. There can be, however, many other 
situations for bosons [7]. Generally, we may take rii bosons to reside in one orbital, (fi{r), 
n2 bosons to reside in a second orbital, </72(r), and so on, distributing the N atoms among 
Uorb > 1 orthogonal orbitals. At the other end to the Gross-Pitaevskii approach lies the 
situation where each boson in the optical lattice resides in a different orbital, i.e., riorb — 
More formally, the multi-orbital mean-field wavefunction for N interacting bosons is the 
following single configuration wavefunction [7] 

^(ri, r2, . . . , riv) = 5 {(pi{ri)(p2{r2) ■ ■ ■ 93iv(rjv)} , (1) 

where S is the symmetrization operator. Note that the Gross-Pitaevskii approach is a 
specific case of Eq. (1) where all orbitals are alike [7,8]. In order to find the ground state of 
the many-bosonic system with the ansatz (1) one has to minimize the energy < ^^lifl^ > 
with respect to ^, Eq. (1). This results in a set of Uorb coupled, non-linear equations for 
the Uorb orbitals that have to be solved self- consistently. To fulfill the variational principle, 
we should search for the energy minimum of H with respect to: (i) the shape of the self- 
consistent orbitals (pi{r), (ii) the occupation of each orbital and (iii) the number Uorb of 
different orbitals. 

The multi-orbital mean- field wavefunction (1) has been successfully employed and led us 
to the prediction of macroscopic fragmentation of repulsive condensates in the ground and 
excited states, see Ref. [8]. In macroscopic fragmentation, a large number of atoms reside 
in a small number of orbitals. Specifically, two or three orbitals were considered in Ref. [8] 



within Eq. (1). What we have found out is that, when the number of orbitals is of the order 
of (number of atoms), the multi-orbital ansatz (1) is physically very relevant for various 
situations of atoms in optical lattices, see below. At the same time, we can now handle 
practically any number of atoms and sites in optical lattices. Specifically here, we employ 
more than 100 bosons/sites, i.e., the number of orbitals is Uorh > 100. More details will be 
reported elsewhere. 

In this work, we concentrate on ID optical lattices, V{x) — Vosm'^{kx), where k is the 
wave vector. As usual, periodic boundary conditions are assumed, thus we employ a super- 
cell of a large number Ny^ of potential wells (unit cells). The strength of the inter-particle 
interaction is expressed via the dimensionless parameter 7 = mg/h'^n, where m is the mass 
of the atoms, n is the density, and g is related to the scattering length and the transverse 
harmonic confinement [10]. Optical potential depths, Vq, and energies will be expressed in 
terms of the recoil energy. 

We begin our study with bosons in a deep optical lattice, Vq — 25Eji, with commensurate 
filling of one atom per lattice site. This familiar SF to MI phase transition is to serve as 
a test tool for our multi-orbital mean-field approach. The MI phase in this system will 
be denoted by MI(1). In practice, we minimized the energy of the Hamiltonian H with 
respect to the ansatz (1). For weakly-interacting bosons we found that all bosons reside in 
one orbital, namely, that the Gross-Pitaevskii approach provides the lowest energy. This 
state, naturally, corresponds to the SF phase. For stronger interactions, the situation for 
the ground state changes completely: each and every atom occupies now a different self- 
consistent, orthogonal orbital! This, as we shall see below, is the MI(1) phase. In other 
words, we employ the ansatz (1) with riorb — N equations to describe the MI(1) phase. In 
Fig. lA, we plot the energy per particle e of the SF and Ml(l) states as a function of the 
dimensionless parameter 7 for N^j — 102 sites. The crossing of the energy curves indicates 
the critical value for which the phase transition from the SF to the MF(1) phase occurs. It 
is found to be 7c = 0.007777(3). It is instructive to translate this value to the "language" of 
the Bose-Hubbard model, i.e. express the corresponding U/ J (on-site interaction divided by 
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hopping) in terms of Vq/Er and 7, see [5,6]. The resulting critical value corresponding to our 
7c is readily found to be {U/ J)c = 3.855(7) — 2 • 1.927(8). It is very encouraging to find out 
that, although our approach is mean-field in real space, it finds a value which is in excellent 
agreement with density matrix renormalization group and other many-body calculations of 
the ID Bose-Hubbard model, see Ref. [9] and references therein. We remind in this context 
that the mean-field critical value of the ID Bose-Hubbard model is {U / J)c ~ 2 • 5.8, rather 
off this value, see [5] and references therein. 

How do the corresponding wavefunctions (orbitals) look hke? In Fig. IB we show the 
SF phase orbital, i.e., the corresponding Gross-Pitaevskii orbital. It, of course, extends 
throughout the optical lattice, representing a coherent state. In Fig. IC shown are a few 
adjacent orbitals describing the MI(1) phase. We remind that these orbitals are obtained 
self- consistently as the solution of a coupled system of Uorh = 102 non-linear equations, 
following ansatz (1). Namely, no preliminary assumptions are made regarding their shape. 
Due to the translational symmetry, a set of equivalent orbitals is obtained, each centered 
around one lattice site. It can be proved that, for vanishing inter-particle interaction, these 
orbitals approachthe lowest-band Wannier functions. This suggests a very appealing physical 
meaning to the ansatz (1). It describes the MI(1) phase with orbitals that can be interpreted 
as boson-dressed Wannier functions. 

With a successful mean-field real-space description of the standard SF to MI transition, 
we set in to exploit more advantages of our method. The next system we considered is the 
SF to MI(1) phase transition in shallow optical lattices (here we took Vq = E^), i.e., for a 
system of strongly-interacting hosom,. For shallow optical lattices the MI(1) orbitals become 
more diffusive, as can be seen in Fig. ID, and extend beyond the next-nearest neighbor 
sites. Accordingly, the corresponding spatial density p becomes flatter - compare Figs. IC 
and ID. Also, notice the negative values of the orbitals in the nearest neighbor wells which 
ensure their orthogonality. For Vq — Er, the value of the SF to MI phase transition has been 
determined to be 7c = 3.4(9). It is instructive to find this mean- field value to agree very 
well with that found by Biichler et al. in the limit Vq ^ by employing the sine-Gordon 
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problem [6]. 

Aiming at exploring more physical situations in optical lattices, we considered another 
scenario of commensurate filling - with two atoms per lattice site. For weak interaction, 
the ground state is, of course, the SF phase where all bosons reside in a single orbital. 
As the interaction increases, a phase transition to the MI phase, which will be denoted 
here by MI(2), occurs. We obtained the MI(2) phase as the ground state of the system by 
ansatz (1). Technically, this is done by determining the Uorb — N/2 self-consistent orbitals 
minimizing the energy with 2 bosons per orbital. As 7 is further increased, it passes another 
critical value, and we find that the ground state of the system underwent a second, new 
phase transition (see Fig. 2A) . Now, there are still two atoms in each site but they reside in 
different orbitals, in contrast to the MF(2) phase where they reside in the same orbital. We 
denote this new MI phase by MF(1,1) (for obvious reasons). In real space, the two atoms 
in each site localize at the borders of the site, leading to minima in the spatial density. 

Figs. 2B and 2C present the wavefunctions (orbitals) of the standard and new MI phases, 
MI(2) and MI(1,1), respectively. The appearance of on-site minima in the corresponding 
spatial density is a clear fingerprint of the new MI phase, differentiating it from the standard 
MI(2) quantum phase. Differences between the standard and the new MI phase should also 
manifest themselves in other observable quantities. For instance, the momentum profile 
of the MI (1,1) phase should be flatter with respect to that of the MI (2) phase. Similarly, 
absorption images resulting from a sudden release of the atoms from the optical lattice in 
the MI(2) phase, where each site houses 2 coherent atoms, should have a slightly sharper 
interference pattern with respect to the MI(1,1), where each site houses 2 atoms in orthogonal 
states. The transition of MI(2) to the quantum phase MI(1,1) when the interaction is 
increased can be interpreted as an on-site transition to the Tonks-Girardeau regime of the 
MF(2) phase. The results obtained so far are straightforwardly extended to more MI phases 
with more atoms per site. For instance, the MI phase with three atoms per lattice site, 
which we denote by MI(3), would eventually end up as the new MI phase, MI(1,1,1), where 
all three atoms reside in different orbitals, and two on-site minima now appear in the spatial 
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density. 

So far, we successfully applied the ansatz (1) to the ground state of atoms in optical 
lattices. This has been obtained by solving a coupled system of norb non-linear equations for 
the orbitals. It is natural to ask whether the ansatz (1) can also provide physical information 
on self- consistent excited states. We remind in this respect that ansatz (1) boils down to the 
Gross-Pitaevskii approach if all orbitals are alike [7,8]. The self-consistent excited states of 
the latter - solitons and vortices - are well-known and have been observed for condensates, 
see, e.g., Ref. [11] and references therein. In optical lattices, the self-consistent excited 
states are the non-linear Bloch bands, the lowest of which is quadratic in k and has zero 
gap [12]. By definition, in the self-consistent excited states of the Gross-Pitaevskii approach 
all bosons reside in the same, higher-energy orbital. In contrast, the flexibility of putting 
atoms in different orbitals suggests that Eq. (1) can also be employed to describe many 
self- consistent excited states. 

In the following, we would like to employ ansatz (1) to understand the nature of the gap 
in the MI phase of strongly-interacting atoms in optical lattices. For this, we recall that 
for weakly-interacting bosons in the MI(1) phase the lowest-in-energy excitation is obtained 
by moving one boson from its site to a neighboring site, see, e.g., Refs. [1,5]. But, what 
happens when the interaction between bosons increases, entering even deeper to the MI 
regime? In that case, the energy of this excited state increases substantially due to the on- 
site interaction between the two bosons. We, therefore, employed ansatz (1) and searched 
for the low-lying self-consistent excited states of the system. We found many such states, 
where, e.g., two atoms reside in a single delocalized orbital, or where two atoms in the 
same site occupy different orbitals locahzed at the borders of this site (similarly to Fig. 2C). 
For sufficiently strong interaction, the excitation lowest in energy does not accommodate 
two bosons in the same site or orbital. Rather, it emerges as an on-site excitation - see 
Fig. 3. It is instructive to compare the energy of the self-consistent on-site excitation with 
the excitation gap of a single well in the harmonic approximation given by \/AVoEr [5] . For 
the parameters used here, 7 = 76.64, Vq — 2bEji and A^^ = 102, we find this value to be 
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0.90-\/4T^£^, slightly lower than the bare-well gap. The difference comes from the strong 
repulsion between atoms which lowers the gap in comparison to the interaction-free problem. 

In conclusion, weakly- to strongly-interacting cold bosonic atoms in deep to shallow 
optical lattices have been described by a multi-orbital mean-field approach in real space, 
giving the wavefunction of cold atoms in the lattice. With it, we described quantum phases 
and phase transitions of cold bosonic atoms in ID optical lattices not accounted for so 
far. As the inter-particle interaction is increased, on-site excitation becomes the excitation 
lowest in energy. The employment of self- consistent mean- field orhitals has been shown 
to provide an accurate value of the SF to MI transition in ID. The findings demonstrate 
the wide potential of our multi-orbital ansatz for cold bosonic atoms and motivate concrete 
studies in higher dimensions. Finally, the predictions obtained are within reach of (effective) 
interaction strengths presently employed in experiments, where values up to 7 = 200 have 
been realized [3] . In addition if needed, Feshbach-resonance techniques can be employed to 
further increase scattering lengths. 
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FIG. 1. Quantum phase transition from the SF to MI(1) phase and corresponding orbitals 93 and 

densities p (shifted lower curves) for weakly-interacting bosons in optical lattice with Vq = 25Eji 

and = 102 sites. Orbitals and densities are normalized (on the segment of length 27r) and are 

plotted against the site index "i" . The optical lattice is illustrated for guidance by the background 

sinusoidal curve. (A) The phase transition is described by the intersection of the SF and MI(1) 

energy per particle curves, e, which occurs at 7c = 0.007777(3). (B) Orbital and density of the 

SF phase for 7 = 0.00776002. (C) Orbitals and density of the MI(1) phase for 7 = 0.00777731, 

slightly above 7c (shown are 3 adjacent orbitals). (D) Orbitals and density of the MI(1) phase for 

strongly-interacting bosons in shallower optical lattice with Vb = Er for 7 = 3.491, slightly above 

the corresponding 7c = 3.4(9) (shown are 3 adjacent orbitals). 
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FIG. 2. Quantum phase transition from the MI(2) to MI(1,1) phase and corresponding or- 

bitals if and densities p (shifted lower curves) for weakly-interacting bosons in optical lattice with 

Vo = 25Eff and Ny^ = 51 sites. Orbitals and densities are normalized (on the segment of length 

27r) and are plotted against the site index "i" . The optical lattice is illustrated for guidance by the 

background sinusoidal curve. (A) The phase transition is described by the intersection of the MI(2) 

and MI(1,1) energy per particle curves, e, which occurs at jc = 12.7(6). (B) Orbitals and density 

of the MI(2) phase for 7 = 12.7418 (shown are 3 adjacent orbitals). (C) Orbitals and density of 

the MI(1,1) phase for 7 = 12.7609, slightly above 7c (shown are 6 adjacent orbitals). Notice the 

on-site minima in the density of the MI(1,1) phase in comparison to the standard MI(2) phase. 
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FIG. 3. On-site self- consistent excited state of the MI(1) phase in optical lattice with Vq = 25£'ij 
and Nyj = 102 sites for strongly-interacting bosons, 7 = 76.64. Orbitals and density (shifted lower 
curve) are normalized (on the segment of length 2tt) and are plotted against the site index "i". 
The optical lattice is illustrated for guidance by the background sinusoidal curve. Shown are the 
excitation site and 4 adjacent orbitals. For the present 7, this is the excitation lowest in energy, 
accommodating the energy gap of the MI(1) phase calculated here to be 0.90V4Vo-E'r. 
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